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Abstract 

We investigate tunneling of excitations across a potential barrier separating two spin-1 Bose- 
Einstein condensates. Using mean-field theory at absolute zero temperature, we determine the 
transmission coefficients of excitations in the saturated magnetization state and unsaturated mag- 
netization states. All excitations, except the quadrupolar spin mode in the saturated magnetiza- 
tion state, show the anomalous tunneling phenomenon characterized as perfect tunneling in the 
low-momentum limit through a potential barrier. The quadrupolar spin mode in the saturated 
magnetization state, whose spectrum is massive, shows total reflection. We discuss properties 
common between excitations showing the anomalous tunneling phenomenon. Excitations showing 
perfect tunneling have a gapless spectrum in the absence of the magnetic field, and their wave 
functions in the low-energy limit are the same as the condensate wave function. 
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I. INTRODUCTION 



Almost ten years ago, an interesting property of an interacting Bose-Einstein condensed 
system was predicted: a potential barrier is fully penetrable for an excitation of a Bose- 
Einstein condensate (BEC) in the low-energy limit Since a potential barrier generally 
inhibits transmission of a particle and leads to total reflection in the low-momentum limit, 



this perfect transmission in the BEC, termed anomalous tunneling [2|, has attracted much 
attention and has been studied extensively and intensively j2-15|. 

The following crucial facts for the anomalous tunneling of the excitation in a BEC (i.e., 
the Bogoliubov excitation) have been discussed: the fact that the wave function of the Bo- 
goliubov excitation has the same form as the condensate wave function in the low-momentum 
limit 16| and the fact that the wave function of its excitation with small but finite momentum 

n 

behaves as the condensate one with a small supercurrent [9|], and so on. 

However, these facts were found through studies of a specific system (i.e., a scalar BEC) 
so that discussions and understandings of anomalous tunneling are restrictive. If we can 
find another system showing anomalous tunneling, it will lead to a deeper understanding 
of anomalous tunneling. An efficient approach in this direction is to examine tunneling 
phenomena in a BEC of particles with spin-1 degrees of freedom (the so called spin-1 spinor 
BEC). 



The spin-1 spinor BEC of dilute atomic^ gases was realized in 1998 16|, ll7|. In the same 



year, Ohmi and Machida [18[ and also Ho [19] determined ground states and excitations in 
the spin-1 spinor BEC through mean-field theory. In this BEC, there are two phases for 
ground states, a ferromagnetic phase and a polar phase, each of which has three-types of 
excitations. For instance, the ferromagnetic phase has the Bogoliubov mode whose energy E 
is linear in the momentum |p|, the transverse spin mode with E oc p^, and the quadrupolar 
spin wave whose energy is massive. 

The following questions will be answered through studies of the spin-1 spinor BEC: Is the 
linear dispersion essential for anomalous tunneling of excitations? (Does an excitation with 
E (X tunnel through the potential barrier without refiection in the limit |p| 0?) Do 
excitations that exhibit anomalous tunneling always have gapless spectra? If we can deduce 
the essence of anomalous tunneling from the answers to these questions, then it will useful 
for predicting tunneling properties in other systems. 
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FIG. 1: (Color online) A schematic picture of the tunneling problem of excitations in a Bose- 
Einstein condensate against a potential barrier. 

In this paper, we thus investigate tunneling problems of excitations in the spin-1 spinor 
BEC. An example of this problem is schematically shown in Fig. [H An excitation runs 
toward a potential barrier from one side and is transmitted through and reflected against the 
barrier. Excitations are scattered off not only the potential barrier but the condensate wave 
function deformed by the potential barrier. The density of the BEC is depleted only around 
the potential barrier because of the nonlinear effect due to the interparticle interaction. We 
study both the ferromagnetic and the polar phases, but we focus rather on the saturated 
magnetization state |m| = 1, and an unsaturated magnetization state < |m| < 1, which 
are, respectively, associated with the ferromagnetic and the polar phases in the absence of the 
magnetic field. Here, |m| is the normalized magnetization, where the magnetization is fully 
polarized when |m| = 1. As well as transmission coefficients, we examine wave functions of 
excitations in the low-energy limit to study a characteristic specific to excitations showing 
anomalous tunneling. 

Table [T] summarizes our main results in this paper. We find that all excitations but the 
quadrupolar spin mode show anomalous tunneling in the low- momentum limit. We also 
find that these excitations have the following two features: (i) the energy gap is zero in the 
absence of external fields and (ii) the wave functions of excitations in the low-momentum 
limit have the same forms as condensate wave functions. 

Using a (5-function potential barrier, we explicitly extract differences between the trans- 
verse spin mode and the quadrupolar spin mode in the saturated magnetization state, both 
of which have spectra with a term proportional to p^. It is worthwhile to compare these two 
modes through an exactly solvable model, since it gives a clue toward understanding how 
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TABLE I: Transmission coefficients T in the low-momentum limit |p| — >• 0. 
Saturated magnetization state (|m| = 1) 



Bogoliubov excitation lim T = 1 

transverse spin excitation lim T = 1 

IpI^o 

quadrupolar spin excitation lim T = 

IpI^o 



Unsaturated magnetization state (0 < |m| < 1) 

Bogoliubov and longitudinal T = 1 

spin excitation (Out-of-Phase) IpI~*'0 

Bogoliubov and longitudinal Wyo. T = 1 
spin excitation (In-Phase) IpI-^'O 



transverse spin excitation lim T = 1 

|p|->0 



perfect tunneling or perfect reflection is induced. A scalar BEC is presumably the simplest 
system in the sense that it has no internal degrees of freedom. However we note that the 
transverse spin wave of the spin-1 spinor BEC is the simplest mode showing anomalous 
tunneling, rather than the Bogoliubov mode in the scalar BEC; the transverse spin wave 
is described by a single-component Schrodinger-type equation while the Bogoliubov mode 
obeys a two-component Schrodinger-type equation. 

Through our studies, the features of anomalous tunneling of the Bogoliubov mode in 
a scalar BEC will be clear. Our results will also lead to further understandings of spin 
excitations in spinor BECs. Moreover, the tunneling problem of excitations is an important 
issue in ferromagnets 20|, l2l| . Beyond the field of cold atoms, the present study leads to a 
further understanding of the spin-wave tunneling in ferromagnets. 

This paper is organized as follows. We review the mean-field theory of the spin-1 BEC 
in Sec. [Tll where the Gross-Pitaevskii-type equation and the phase diagram of the ground 
state are derived. The Wronskian and the energy flux are also discussed to determine the 
transmission coefficients. Section IIIII presents the tunneling properties of excitations in the 
saturated magnetization state. Section HV] presents the tunneling properties of excitations 
in an unsaturated magnetization state. In Sec. |Vl we investigate the tunneling problems 
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of excitations across the (5-function potential barrier in the saturated magnetization state. 
Using an exact solution of the wave function, we derive the transmission coefficient of the 
transverse spin mode. We also analytically discuss perfect reflection of the quadrupolar 
spin mode using the 5-function potential barrier. In Sec. Wl\ we examine the dependence of 
the tunneling properties on the spin-dependent interaction strength and the linear Zeeman 
effect. Appendix lA] serves as a summary of the eigenmodes and Appendix [B] presents the 
tunneling properties of the unsaturated magnetization state in the integrable case. 



II. MEAN-FIELD THEORY OF SPIN-1 BEC 



A. Gross-Pitaevskii-type equation and ground state 



We summarize the mean- field theory of a weakly interacting spin-1 Bose system developed 
by Ohmi and Machida 18[ and Ho 19| in order to fix notations. We start with a Hamiltonian 
of a spin-1 Bose system described by 



H 



i 

+ 1/ ^r5^^1(r)^Kr)^,(r)^,(r 



+ 



Cl 



(1) 



where the sums with respect to i, j, i' and j' are taken over the magnetic sublevels ±1 and 
0. ^i(r) is the Bose field operator satisfying [^i(r), ^](r')] = 5ij5{r-v'), [^i(r), ^j(r')] = 0, 
and [^^|(r),^](r')] = 0. The sum with respect to k is taken over the Cartesian coordinates 
X, ?/, and z\ m is the mass, g is the Lande's g factor, /xb is the Bohr magneton, and B{> 0) 
is the strength of a magnetic field (with the z axis being taken as the direction of the 
magnetic field). The interaction strengths are given by Cq = 47r/i^(2a2 + Oo)/(3^«), and 
Cl = 47r/i^(a2 — ao)/(3m), where ast^t is the s-wave scattering length for a collision process 
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with total spin S'tot- The spin matrices of the spin-1 system are given by 



V2 



V2 



0-10 
1 -1 
1 



1 

0-1 



(2) 



1 1 
yO 1 Oy 

The field operator obeys the equation of motion ihd'^j/dt = [^j, H]. Within the mean- 
field theory (^j) = $^6"*'^^*/'', the Gross-Pitaevskii(GP)-type equation of the spin-1 BEG 
for stationary states is given by 

(3) 
(4) 

where n = -h^V'^/{2m) + - n + cqh and F+ = F*_ = y/2Re\{^\^^Q + We 
define also n = + |$oP + I'^'-iP and = IJLjHiP - I'^-iP- We use fij = /i for j = ±1 

and 0, satisfying the condition 2fio = /x+i + [22]. 

In the following, we consider the phase diagram of the ground state in spatially uniform 
systems in the presence of magnetic fields. It is easy to see that a fully polarized state 



n T gf^BBj $±1 + c,{F^%/V2 ± F,$±i) = 0, 
n% + Ci(F+$+i + F_$_i)/v^ = 0, 



($+i,$o,$-ir = (v^,0,0) 



(5) 



is a solution to (|3]) and (jl]). The total energy per unit volume S = —fin + cqu^ /2 + £' , with 



(6) 



for ([5j) is given by S' = c\r? j2 — q^bB. A state may be characterized by the magnetization 
M = (M^, My, M') with M'' = Ej=±i,o "^i^ij^i k = x,y, z. When we denote by m 
the normalized magnetization [i.e., m = M/(/iB^)], |m| = 1 follows for ([5]). We name the 
state with |m| = 1 the saturated magnetization state. Another solution to ([3]) and (jl]) 
is given by 



($+i,$o,$_i)^ = (0,y^,0)T, 
which has S' = and m = 0. When (0 <)g^BB < cin, the state 

(<l>+i, $0, = ^/^i±VlTA, 0, v^r^)^. 



(7) 



(8) 
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with A = g^BB/ciu is the third solution to ([3]) and (jl]). This solution has £' = 
— ((7/iB-B)^/(2cin) and |m| = A. We name a state with |m| G [0, 1) an unsaturated magneti- 
zation state. Among (j7]) and ([HD, the solution with the lowest £' yields the ground state. 
The ground state is given by the saturated magnetization state when cin < q^bB, while 
the state is realized as the ground state when (0 <)gfiBB < Ciu. 

When B = and ciK 0, the saturated magnetization state having |m| = 1 is called the 



ferromagnetic phase 19|. When B = and ci > 0, however, the unsaturated magnetization 
state having |m| = is called the polar phase |l9j]. Note that the degeneracy of the ground 
state relating to the gauge symmetry and the spin rotational symmetry exists when B = 
0. A configuration of an order parameter (the condensate wave function) is transformed 
into another through the gauge transformation e*^° and the spin rotations U{a,P,'j) = 
^-■iS^'a^-isyp^-iS''^ _ ferromagnetic phase, the order parameter is given by [l9(] 

(<l>+l,<l>0,$-l)^ 

=e^^°W(a,/3,7)(v^,0,0)^ 

=v^e^(^«-^) (e-^"cos2^,y2cos^sin^,e*"sin2^y . (9) 
In the polar phase, the order parameter is given by [loj ] 

($+l,$0,$-l)^ 

=e*''«W(a,/3,7)(0,v^,0)T 

=yne^^o (^--^e-*"sin/3,cos/3,--^e^"sin/3^ . (10) 



When one takes the z axis as the spin-quantization axis [19|, a configuration of the or- 
der parameter in the polar phase is reduced to ($+i, $o, ^-i)"*" = (0,v^5 0)"'". The 
spin quantization axis for ($+1, $0) '^'-i)^ = (— a/?t./2,0,a/?t-/2)'^ is the x axis, and for 



($+1, $07 '^'-i)^ = (a/7^/2, 0, A/n/2)'^, it is the y axis. As a result, these configurations are 
equivalent, and the state ([7]) is equivalent to the state (j8]) when 5 = 0. 

B. Equation of excitations and transmission coefficient 

In order to study small oscillations of the system around equilibrium, we consider a 
small fluctuation 0i(r, t) deviated from the condensate wave function $i(r), where 0i(r, t) is 
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regarded as the classical field (c-number). The equation of excitations is given by 



ih- 



-"±1 



dt 



^±1 



+ ^£Ll0Tl + (^+ll-l+Cl'^o)0Tl' 



ih- 



dt 



'H + cin + ) 00 + (co$g + 2ci$+i<l>_i: 



+ J2 iP*J^ + R^ 



■Q,rar'm)i 



(11) 



(12) 



m=itl 



where P±i = (cq + Ci)<l>S<f ±i + 2ci^*^,%, = (cq ± ci)$,%, and Pjf = (cq ± ci)$*%. 

A spatially constant quantity is important to determine the transmission and reflection 
coefficients of excitations. We derive this quantity from Eqs. f|TT]) and f|T2|) . By assuming 



v\{r))*e''^^^l^], with a label ^ 



the Bogoliubov-type wave function (l)i{r,t) = ^[M'(r)e~*^'*/^ 

I 

of an eigenmode, we can obtain the Bogoliubov-type equation for the spin-1 BEC from Eqs. 



Till and (1121) . With the use of this equation, we have the following relation: 



i=0,±l 



<4'(r)(,4(r))--t.f(r)(!.'(r))-] 



.__v.w,,w, 

where the Wronskian Wi ;/(r) is given by 

W,,(r)= J2 [(^r(r))*Vn^(r)-«Kr)V(«r(r))' 

i=0,±l 



(13) 



(14) 



Here / and /' denote labels of eigenmodes. When eigenmodes / and /' have the same energy 
E = = E^\ one finds Wi i/(r) = const. Hereafter, we discuss the Wronskian W;^; (i.e., 
/' = /) to determine the transmission and reflection coefficients. 

This Wronskian of the Bogoliubov equation is related to the time-averaged energy flux 
Q]. An equation of the total energy 8 per a unit volume is given by 

where the energy flux Q is defined as 

-- E 

2m dt 



Q 



(15) 



(16) 
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Let us introduce a time-averaged quantity (X) given by 

{X) = - [ Xdt, (17) 
^ Jo 

where r is the time given by r = 2TThn/E {n is integer) and E is the excitation energy. 

We consider the time averaged energy flux of the eigenmode /, assuming that the small 
fluctuation deviated from the order parameter $j(r) is given by the Bogoliubov form, i.e., 

= + wje-'^'*/^' - {v^^Ye'^''^^ (18) 

As a result, the time averaged energy flux (Q^) of the eigenmode / is given by 

W = - ^W,,,. (19) 

Since 'Wi^i is spatially constant, the time-averaged energy flux of an eigenmode / is also 
independent of the position. 

We give the transmission and reflection coefficients using this spatially constant quantity. 
In this paper, we study one- dimensional tunneling problems, assuming the normal incidence 
of excitations against a potential wall. We name the incident side of the tunneling problem 
the left side and label it as L. We label the right side as R. The spin-1 spinor BEC has some 
types of excitation because of the internal degrees of freedom. In particular, as mentioned 



below (and also in Appendix |A]), three modes in the uniform system are well known 18l. Il9|. 
which may characterize the incident, reflected, and transmitted waves. Let us suppose that 
the incident mode is one of these three excitations in the spin-1 BEC, which we name a = I. 

n 

We name the other two modes cr = II and III [23]. A wave function far from the potential 
barrier on the left side is described by a sum of the incident and reflected waves: 

«i,L(r), i;+i,L(r), <L(r), <L(r), «-i,L(r), t'ii,L(r))'' 

o-=II,III 

The wave function far from the potential barrier on the right side is described by the trans- 
mitted wave: 

«i,R(r), <i,R(r), <R(r), v^^ir), ui.^ ^(r), ^(r))^ 

~^ 5^ ^^6*^°^^ X (''^+l,R5 ^+l,R5 ^0,R5 ^0,R' ^-l,R' ''^-I.r)'^- (^1) 

cr=I,II,III 
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j=0,±l 






{Qa,L) 




E (i<Lr+i<Lp 

i=0,±l 


(for (7 


= 11 and III), 


{Qa,R) 


2m 




% (for a 


= 1,11 and III) 



Here, and i^, respectively, are the amplitude reflection and transmission coefficients of a 
mode (7. The time-averaged energy flux of a mode a in each asymptotic regime is given by 

(22) 
(23) 
(24) 

The transmission and reflection coefficients To- and of a mode a are, respectively, 
defined by the ratios of transmitted energy flux to incident energy flux, and also the ratios 
of reflected energy flux to incident energy flux. Using this deflnition, we obtain and 
as 

^.,R^_E^(I<RP + I<RP) 

^'^^^l'e^KlP+KlP) ^^^^ 

i=0,±l 

,0- \2\ 



i=0,±l 

""^I,L.E (KlP+KlP)''^' ■ 



^.=i-^^^^W^^::r^k^r- (26) 



i=0,±l 

According to (Qi^i) = {Qi,r), one has Ea(^<^ + -^<t) = 1- 

These are a general framework for the tunneling problem of excitations in the spin-1 
spinor BEC. In the next section, we examine the tunneling probability of excitations in the 
spin-1 BEC focusing on the saturated magnetization state. After that, we investigate the 
unsaturated magnetization state. 



III. SATURATED MAGNETIZATION STATE 



The saturated magnetization state is realized when cin < q^bB. Its condensate wave 
function of the uniform system is given by ($+1, $o; $-1)^ = (\/^, 0,0)^. Substituting 
the above wave function into the GP-type equation, one obtains the chemical potential 
/X = c+n — q/ibB, where c+ = Cq + Ci. It is expected that wave functions $0 and $_i do 
not appear near a potential barrier, since a modulation of condensate wave functions costs 
kinetic energy; thus, we assume the following condensate wave function near the potential 
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FIG. 2: Total density profiles of the saturated magnetization state in the presence of the rect- 
angular potential barrier (j29p . with the various potential heights V\y = Vb/(CfC+'^) and with 
Vi^ = Vk = 0. 



wall: 



(<l>+i,<l>o,$-ir = (v^0(x),O,O)^ 



(27) 



We have checked the validity of this assumption through the numerical calculations of the 
GP-type equations ([3]) and (jl]). Substituting fl271) into the GP-type equation, we have 







+ Kxt - c+n (j) + c+n(j)-^. 



(28) 



2m dx^ 

We notice that this equation is the same as the GP equation of the scalar BEG. Figure [2] 
shows the density profiles for the rectangular potential barrier given by 



V{x) 



for X < 0, 

for < X < 6/V2, 

for ^f/A/2 < X. 



(29) 



The healing length of the ferromagnetic state is given by .^f = h/ y/mc+n. The magnetization 
is proportional to this profile. As in the schematic of Fig. [H the condensate density is 
depleted around the potential barrier. 

Using this wave function for the condensate, we obtain the following equations of fluctu- 
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at ions: 



ih 



dt 



ih- 
in- 



^0 



dt 



2m dx^ 
2m dx"^ 

d^ 



+ Kxt + c+n{(j) - 1) 0+1 + c+n(f 
+ Kxt + c+?2(0^ - 1) + g/J^BB 00, 



1 + 



(30) 
(31) 
(32) 



dt V 2m dx"^ 
where c_ = cq — Ci. 

We use the Bogohubov transformation </)+i = Uj^iexp {—iEt/h) — v*^iexp (iEt/h) and 
introduce two functions, S^i = u+i + v+i and G+i = m+i — ti+i, for the magnetic sublevel 
+1. To obtain the transmission coefficient for the spin component +1, we solve the following 
equation: 

d^ 



EG 



+1 



ES+, 



2m dx"^ 



+ Kxt(a;)+c+n(0'-l) 



+ V'cxt(a;) + c+n 



2m dx"^ 

imposing the boundary condition 

I = [exp (ikx) + r exp (—ikx)] i | +aexp(/tx) 



a+i 



S+i 



-texp{ikx) I I +6exp(— Kx) j ^^'^ 



(33) 
(34) 

for a; — )■ — oo, (35) 
forx — > +CXD, (36) 



where k and k are, respectively, given by hk = J2m[^{c+ny + E^ — c+n] and = 



2m\^ (c+n)^ + i?^ + c+n]. The excitation of the +1 component becomes the ordinary 
Bogoliubov excitation with E = \/e{e + 2c^n), where e = h'^k'^ /(2m), as noted in Ap- 
pendix |X1 a+i and /3+i are also given in Appendix |X1 
For excitations 0o and we solve the equations 

d^ 



2m dx'^ 

d^ 



+ Kxt + c+n( 



2m dx"^ 

imposing the boundary conditions 

(f)j — > exp (ikx) + r exp {—ikx) 
— T-t exp (ikx) 



+ Kxt - c+n + c_n0 + 2gfiBB 



for x — )■ — oo, 
for X — )■ +00, 



(37) 
(38) 

(39) 
(40) 
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FIG. 3: (a) Transmission coefficients T and (b) corresponding pliase sfiifts 9t (defined by tlie 
argument of the amplitude transmission coefficient t) of tlie Bogoliubov mode with the spin com- 
ponent +1 in the saturated magnetization state through the rectangular barrier (j29p . We used 
^ = = 0, and define V\y = V\^/{c+n). We use parameters oq : 02 = HO : 107 following 
Ref. Q- 

with j = and —1. Here the momentum hk for the spin component is given by hk = 
yJlmiE — gj2-BB) and for the spin component —1, by hk = A/2m [E + 2(cin — gfi-^B)]. Ex- 
citations of the magnetic sub levels and —1 are, respectively, associated with the transverse 
spin mode and the quadrupolar spin mode with E = e + q^-bB and E = e — 2cin + 2gfiBB, 
as noted in Appendix |X1 Only the quadrupolar spin excitation has the non-zero inherent 
energy gap Ab=o = 2|ci|n for A; — )■ when B = 0. We use the term "the inherent energy 
gap" as the energy gap in the case when B = 0. The other excitations in the saturated 
magnetization state have no inherent energy gap (i.e., A^^o = 0). 

We show the momentum dependence of the transmission coefficients and phase shifts 
for all three excitations in Figs. [3] and |H We find that the Bogoliubov mode and the 
transverse spin wave show the anomalous tunneling phenomenon (i.e., perfect tunneling 
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FIG. 4: (al) Transmission coefficients T and (bl) corresponding phase shifts 9^ of the transverse 
spin wave mode in the saturated magnetization state. (a2) Transmission coefficients T and (b2) 
corresponding phase shifts 6t of the quadrupolar spin wave mode in the saturated magnetization 
state. The barrier (j29l) and the parameters oq and 02 are the same as those in Fig. [3l 

in the long- wavelength limit), although these two excitations have different spectra. Phase 
shifts, defined as the argument of the amplitude transmission coefficient t, of both excitations 
reach zero in the low-momentum limit. Unlike these two excitations, the quadrupolar spin 
mode does not show the anomalous tunneling phenomenon. Its phase shift does not reach 
zero in the low-momentum limit; it reaches tt /2 in this case. 

Interesting characteristics can be seen in the wave functions of the excitations. Figure [5] 
plots wave functions of the Bogoliubov mode in the presence of the rectangular potential 
barrier fl29l) . One finds that the wave function E / {2cj^n) S+i{x) in the long wavelength 
regime coincides with the condensate wave function (i.e., S+i oc $+1). The existence of 
this type of solution in the limit E can be easily confirmed by comparing the GP 
equation ( 128|) with ( 133|) . From Fig. [5l however, one confirms that G^i{x) is absent in the long 
wavelength regime (i.e., limfc^o G^+i = 0). As a result, one finds that (m+i, f+i) oc ($+1, $+1) 
follows. Through studies of anomalous tunneling in the scalar BEG, it was pointed out 
that total transmission in the long wavelength limit is deeply related to a correspondence 
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5S 
1^ 



5S 



fc^ = 0.003 

fc4 = 0.2 

(al) real part 


/ 

(a2) imaginary part 






(bl)real part 


(b2) imaginary part 



-6 - 



FIG. 5: (Color online) Wave functions of the spin component +1 in the saturated magnetization 



state. S{x) = S'+i(j;) Y'^7(2c+n) and G{x) = G+i {x)y^E/ (2c+n) . We use the barrier (p9]) . where 
^ = ^ = and Vb = 3c+n. The parameters oq and 02 are also the same as in Fig. [3l 

between a wave function of the Bogoliubov excitation in the low energy limit and that of 
the condensate |6|]. This characteristic for the Bogoliubov mode in the scalar EEC is the 
same as that for the Bogoliubov mode in the spin-1 EEC 

Let us focus on the other two modes in the saturated magnetization state: the transverse 
spin-wave mode associated with 0o, and the quadrupolar spin mode associated with 
Figure [6] plots the wave functions of these modes. From panels (al) and (a2), we see that 
the wave function 0o in the low-energy limit has the same form as the condensate wave 
function $+1. Comparing the GP equation (l28l) with (l37j) for the limit — )■ 0, we can 
confirm the existence of a solution being proportional to the condensate wave function $+1 
(i.e., 00 oc $+1). As for the spin component —1, there are no solutions proportional to <l>_^.i 
for /c — > 0. From panels (bl) and (b2), we find that the amplitude of 0_i goes to zero as 
the wavelength becomes longer; this leads to the total reflection (i.e., t = and r = — 1) 
in the limit — > 0. This is a characteristic feature of the quadrupolar spin mode, which 
was noted in our earlier work [l3|. We will revisit the quadrupolar excitation in Sec. IVj and 
show that the potential barrier makes the amplitude of its wave function absent in the long 
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FIG. 6: (Color online) Spatial profiles of the wave function, (al) and (a2) plot cpQ] (bl) and (b2) 
plot (p-i- The potential barrier and parameters are the same as those in Fig. [3 

wavelength limit. 



IV. UNSATURATED MAGNETIZATION STATE 

We examine the transmission coefficients of excitations in the unsaturated magnetization 
state, which is realized when (0 <)g^BB < cin. The wave function of its ground state in the 
uniform system is given by Eq. ([8]). Substituting this configuration into the GP equation, 
the chemical potential is given by /x = Cqu. $o is expected to be absent everywhere, for 
the same reason mentioned in the previous section, and hence we assume the following 
configuration: 



($+i,$o,$-i)^ = ($+1,0, $_i 



(41) 



We have checked its validity through the numerical calculations of the GP-type equation 
([3]). Substituting the above configuration into the GP-type equation, we obtain the following 
equation of $-1-1: 







2m dx"^ 



(42) 
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In Fig. [TJ we show the condensate wave functions, number densities of particles, and the 
magnetization for various magnetic fields in the presence of the potential barrier. It is clearly 
seen that the density profile does not depend on the uniform magnetic field, while the profile 
of the magnetization is easily changed by the uniform magnetic field. 

On the basis of the configuration (1411) . we obtain equations of fluctuations, which are 
given by 



ih- 



■>±1 



ih- 



dt 



2m dx"^ 



60 + 2ci$+i$_i( 



(43) 
(44) 



When we introduce the Bogoliubov-type transformation 0j = UjC «'E^*/^_^*g+*'Et/?i ^^d define 



Sj = Uj + Vj and Gj = Uj — Vj, we obtain the following equations: 

ES+i = (n+ + 2c+<I>^i) G+i + 2c_<l>+i$_iG_i, 

EG— I = H—S—i, 



where i-L± 



(45) 
(46) 

(47) 
(48) 



— + \4xt — Cqu =f QfisB + c+$L + c_$|i. In order to determine the 

2m dx^ ^ 

transmission coefficients, we impose the following boundary condition: 



forx — )■ —00, 

























+ 1 

j = 












\G-J 
































G+i 


i=± 















(49) 



forx — )■ +00, 



\G-iJ 



(50) 
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FIG. 7: (Color online) Condensate wave function (a) $+i, (b) (c) the local density and (d) 

the magnetization for various magnetic fields in the presence of the 5-function potential barrier 
Vext{x) = Vh6{x) for Vb/^p = Sam, where = h/ ^mcon. B = gfiBB /{cin). We use parameters 



oo : 02 = 46 : 52 following Ref. 
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where /c(±) and K(±) are given by hk(^±) 



m 



v/(c+n ± C)2 + 4^2 - (c+n ± C) 



and 



(c+n ± C) , respectively. Details of excitations for 
the uniform system [such as the wave numbers of the oscillating and evanescent waves, k± 
and K±, and the amplitudes of the wave function, q;(±), /3(±), o:{±), and /3(±)], are summa- 
rized in Appendix IXI The mode specified by /;;(+) (k(+)) or A;(_) denotes the out-of- 
phase mode or in-phase mode between the spin components ±1. Its energy E(^±^ is given 
by -E'(±) = \^/e(e + c+n =fc C), where e = !fk'^j_^/{2m) and C = a/ (c_r7,)2 + AcQCirfiA^ with 



A = gfi-gB /{ciTi). In Ref. [18|, these excitations are regarded as the collective modes coupled 
with the number density and the longitudinal spin density. These excitations are gapless, 
and the inherent energy gap, A^^o = 0, is zero. These excitations in the low-energy limit 
are regarded as phase modes of condensates, as seen in the Bogoliubov mode, so that we 
name the mode with [-E'(-)]; "Bogoliubov and longitudinal spin excitation [put-of-phase 
(in-phase)]" as in Table [H 

The transmission coefficient of the excitation (po is, however, obtained by solving the 
following equations: 





(uq] 




H 




[Vol 





-2ci<l>+i<l>_i 1 


II 


(uo\ 


1 —'Ho , 







(51) 



where 'Hq 



^2 rf2 



2m dx'^ 

transformation given by 0o 
conditions: 



+ Vcxt — cqu + c+($2_^ -|- $^i), and we used the Bogoliubov-type 



-iEt/h 



WqC"^*'^*/''. We impose the following boundary 




+ re 

















for X 



-oo, 



where k and k are given by hk 




foTX — )• -|-oo. 



(52) 
(53) 



^J2m A/(cin) 



;i - A2) + ^2 - cin 



and hn 



2m 



a/ (cin)2(l — ^2) + £'2 + cin , respectively. (Details of this excitation for the uni- 
form regime are also given in Appendix El ) This mode is associated with the transverse 
spin wave mode with the energy E = >/£(£ + 2cin) -|- {q^bBY- This mode has the energy 
gap qii-qB, but the inherent energy gap is zero (i.e., Ab=o = 0). 
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Figure [8] plots the results for the incident mode being the out-of-phase mode (left panels) 
and the in-phase mode (right panels); Fig. [9] plots the results of the transverse spin wave 
mode. Here the wave number k is scaled as fc^p, where = h/^/mc^. We also find 
that all excitations tunnel through the potential barrier in the long wavelength limit, where 
phase shifts of all modes reach zero. We find that the sum of T and R for the out-of-phase 
(in-phase mode) is unity, as seen in Fig. |8l As a result, when the incident mode is the out-of- 
phase (in-phase) mode, the transmitted and reflected waves are the out-of-phase (in-phase) 
mode itself and do not include the opposite type of mode [i.e., the in-phase (out-of-phase) 



mode] [24 1 . 



As in Sec. IIIIl we saw that the wave function of an excitation has the same form as the con- 
densate wave function when the anomalous tunneling phenomenon occurs in the saturated 
magnetization state. Here, let us examine the low-momentum properties of the excitations 
in the unsaturated magnetization state. Figures [10] and [TTl respectively, show wave functions 
of the out-of-phase mode and the in-phase mode. As for the out-of-phase mode, one finds 
that 5*+! and S^i are, respectively, proportional to $+1 and <l>_i; but G±i is absent in the 
low- momentum limit. Consequently, from Fig. [101 the solution in the long- wavelength limit 
for the out-of-phase mode is given by (u+i, t^+i, u„i, f _i) oc ($+1, $+1, $_i). Because 
of the opposite sign between $-1-1 as in Fig. [TJ this excitation can be confirmed to be the 
out-of-phase mode. 

As for the in-phase mode, one first finds that G±i is absent for the small momentum 
from Fig. [HI Also from Fig. [TTJ however, one finds that Im[S'-i-i] is not absent, although 
Re[5'±i] is proportional to $±1. This is due to the momentum fc^f = 0.001 being not as small 
as the low momentum limit for the in-phase mode. In fact, Im[S'-|-i] has a gradient in the 
asymptotic regime. In the low momentum limit, the gradient of the wave function for the 
tunneling problem should be absent far from the potential barrier, because the boundary 
condition of this problem is given by exp {ikx) + rexp {—ikx) for x <^ — ^p, and texp (ikx) 
for X ^ .^p, so that this condition leads to 1 + r(= const.) for x ^ — .^p and t{= const.) for 
X ^ ^p in the low-momentum limit. To examine the wave function S±i of this mode in the 
low-momentum limit, we return to Eqs. fj45|) and fj47|) . In the limit E ^ (i.e., k — )■ 0), 
Eqs. (gSD and (gTD lead to 

=n±S±i. (54) 
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If we compare ( H2|) with (15^ . a set of linearly independent solutions to (15^ is found to 
be given by $±i(x) and $±i(a;) dx'^^i{x'). The solution 5*^1 is now given by the linear 
combination of these two functions. However, since = =fV1 ± A(= const.) for |a;| ^ ^p, 
the second solution is proportional to x for large This means that the second solution 
does not satisfy the boundary condition of the tunneling problem in the low- momentum limit, 
and is found to be absent. As a result, S±i oc $±i follows. Consequently, the solution in the 
long-wavelength limit is given by f+i, f _i) oc (— — where the 
sign is given by the boundary condition. 

We finally discuss the transverse spin excitation. Figure [12] plots the spatial variation of 
wave functions uq and vq- From this figure, we find that the wave function of the excita- 
tion in the low- momentum limit has the same form as the condensate wave function [i.e., 
— ^/A{uq,Vo) = (<l>+i,$_i)]. This property can be confirmed from flSTl) in the limit /c — )■ 
(i.e., E — q^-qB). In this limit, one finds that Uq and Vq have solutions Mq oc $+i and 
fo oc $-1 by comparing (IST]) with ( l42l) . From those results, (Mo,fo) (*l'+i, "l*-!) follows in 
the low momentum limit. 

In summary, recalling the case of the saturated magnetization state, we found that all 
excitations showing total transmission in the long-wavelength limit have wave functions 
equal to the condensate wave functions. As for the Bogoliubov mode in the scalar BEC, 



Ref. 



25| pointed out the existence of the solution equal to the condensate wave function 



for the zero mode. The zero modes in the spin-1 BEC also have the same property; their 
wave functions in the low momentum limit are equal to the condensate wave functions for 
the spin-1 BEC. 



V. EXCITATIONS IN THE SATURATED MAGNETIZATION STATE IN THE 
PRESENCE OF THE 5-FUNCTION POTENTIAL BARRIER 

In this section, we examine solutions in the presence of the (5-function potential bar- 
rier, focusing on spin modes in the saturated magnetization state. As for the Bogoliubov 
excitation, the property in the current free state was studied in Ref. jl|. 

First, we derive an analytic form of the transmission coefficient for 0o- We assume the 
potential energy Vey±{x) as Vl for x < 0, Vr for x > 0, and Vb5(x) for x = 0. As shown in 
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Eq. (128!) . the GP-type equation in this saturated magnetization state is given by 



2m dx'^ 



$+1 = 0. 



(55) 



Let us introduce the heahng length ^l(r)5 given by ^l(r) = y "^(/U — VL{R))//i. In the asymp- 
totic regime, we have the following forms: $+i = A/(/i — Vl)/c_|_(= y/rvi) for x <^ —^L: 



and $+1 = A/(/i — Vr,)/c+(= ^/n^) for x ^ ,^r,. In the presence of the (5-function potential 
barrier, the condensate wave function is given by 



$+1 



(56) 



/^L tanh [(— X + Xo,l)/^l] for x < 0, 
/riRtanh [(x + a;o,R)/^R] for x > 0. 

Parameters Xo,l and Xo,r are determined by sewing the wave functions at the barrier [i.e., 
<l>+i(+0) = $+i(-0) and <I>Vi(+0) - $'+i(-0) = 2mV^<!>+i{0)/h' ]. 
In this case, the solution to (!37|) is given by 

for X < 0, (57) 
for X > 0, (58) 



^0(3;) =fh{k, x) exp {+ikx) + rfi^{—k, x) exp {—ikx) 
t>o{x) =tfR{k, x) exp {+ikx) 



where fi^{k,x) = — zA;^,ySgn(x) +tanh [(| 

,u)/^v]- Sewing the different pieces of the wave 
functions together at x = [i.e., 0o(+O) = 0o(-O) and 0o(+O) - (j)'o{-0) = 2mVb0o(O)//l^], 
we have the amplitude transmission and reflection coefficients t and r, given by 



t 



\riL-ikiiMil + il-2ikr^^i^iiy 
.'2^kr^lil.il-^k{a + il)-Vl.{a-eK) 



(59) 
(60) 



(r/L-^A;eL)(e^ + ei-2^fcryLeLei) ' 
where ?7l = tanh (xo,l/^l)- In the long-wavelength limit, we have the transmission coefficient 
T = |tp in terms of the magnetization My ^a^riy, which is given by T = 4MlMr/(Ml + 



Mr) 2. This result recovers the result in Ref. 



15| . We conclude that the transverse spin- wave 



mode in the saturated magnetization state shows partial transmission in the long-wavelength 
limit when Ml 7^ Mr. 

In the case where Ml = Mr, we have 

vil + k'e) 



iv - ^^0(1 - ikr]C)'' 

ikiirf - 1) 
{V - ikC}{l - ikr]iy 



(61) 
(62) 
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where t] = rji^ and ^ = = ^r- As a result, total transmission in the limit /c — )■ occurs. In 
the high-barrier limit, we have rj ~ h"^/ [m^Vh), so that the amplitude transmission coefficient 
in the long- wavelength regime is given by t ~ l-l-i/cm^^Vb/^^. In this limit, the transmission 
coefficient T is given by T = |tp = 1 + 0{k'^). We also find (t)Q{x) = up to an overall 

normalization factor in the limit A; — )■ 0. The correspondence between wave functions holds, 
as discussed in Sec. IIIII 

We note that it is easy to extend this one- dimensional tunnelingproblem to the problem 
of the reflection and the refraction of the excitation as studied in [7] . [The barrier has only 
the X dependence (i.e., the condensate wave function also has only the x dependence], and 
the incoming mode has the incident angle deflned by that between the x-direction and the 
incident momentum.) Even if Vl 7^ Vr, the energy spectrum E = /i^fc^/(2m) -|- q^bB is 
independent of the potential barrier. Because of the translational invariance, the momenta 
along the wall on both sides are equal. As a result, the equation of 0o is obtained when we 
replace E in Eq. ( l37|) with fi^k"^/ {2m) + q^bB, where hk^ is the momentum along the x-axis. 
Consequently, one flnds no refraction, and the fc^-dependence of the amplitude transmission 
and reflection coefficients is obtained if one replaces k in ( l59l) and ( I6OI) with k^. 

The 5-function potential barrier is also useful to study the total reflection of the quadrupo- 
lar spin mode in the saturated magnetization state. In what follows, we investigate a sim- 
ple model: the junction of BECs with equal densities separated by the potential barrier 
Vext{x) = Vb6{x). We compare the properties of wave functions 0o with 0_i in the presence 
of this potential barrier. The simple problem with the 5-function potential barrier has re- 

nn 

vealed essential phenomena for the tunneling problem of the Bogoliubov excitation |3[ . 
The implication of the tunneling problem against the 5-function potential barrier is expected 
to be common in problems for any potential barriers with an arbitrary shape. 
The GP-type equation (l28l) normalized by c+n is given by 

= 0, (63) 

where we replaced x/^f with x, and Vt,/ (^c+n) with Vb- In this case, the solution is given by 
0(x) = tanh(|x| + xq), where the parameter xq is determined hj r] = tanh(xo) = — (Vb/2) -|- 
a/ (Vb/2)2 + 1 through the boundary condition at x = [i.e., 0(+O) = </>(— 0) = rj and 
0'(+O) = 0'(-O) + 2l^b0(O)]. From Eqs. ^ and (ES]), an equation of 0^ for j = or -1 at 
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Ixl > can be reduced to 



'2 



Id 



y0., (64) 



where Aq = 1 and A^i = 1 — 2ci with ci being ci/c+. 

Using the fact d(f){x)/dx = sgn(x)[l — 0^(x)], we transform Eq. into the following 

form: 



1- 



A general solution to (!65|) is given by 

0, = «J(fc)i^;(0) + /3J(A;)g;^(0), (66) 

where P^*^ and Q^^, are associated Legendre functions, and a parameter (^j is given by Q = 
(-1 + ^1 + 8Aj)/2 [i.e., Co = 1 and C-i = (-1 + V9 - 16ci)/2]. The index v in coefficients 
and /3J represents L for a; < and R for a; > 0. 

An associated Legendre function Q^^'^°(0) diverges at |x| — > oo because of <^i\x\ — )■ oo) = 
1, and hence we have a condition fi^^{k — > 0) = 0. Sewing the different pieces of the wave 
functions together at x = [i.e., 0j(x = +0) = 0j(x = —0) and 0^(+O)— 0j(— 0) = 2Vb0i(O)], 
we have two conditions: = q;^(A; — )■ 0) = o^{k — )■ 0), and 

ik f i\ ryik I 



d_P^_P^) 
'ct>^r]k^O \ d(j) (f) 



= aj lim lim^ | — ^-^j— ] . (67) 



Here we used the boundary condition of the condensate wave function at x = 0. 

As for the mode j = 0, P^^^^-^ (</>) = (f) holds. As a result, one finds that oq is not generally 
zero from Eq. (l67j) . so that (pQ cc (p follows. We here consider the boundary condition of the 
tunneling problem [i.e., exp {ikx)+r{k) exp {—ikx) for x ^ — 1 and t{k) exp (ikx) for x ^ 1]. 
In the long wavelength limit — )■ 0, we have 1 + r{k — )■ 0) = t{k — )■ 0) = </)(|x| — t- oo) = 1, 
where we set = 1. From this condition, we have total transmission t = 1 and r = 0. We 
also obtain the phase shift of the amplitude transmission coefficient, 6t = —i\n{t/\t\) = 0. 

In the case of the quadrupolar spin mode for -P^^^^]^)(0) 7^ holds. We conclude 
that a_i should always be zero in order to satisfy the condition ( 1671) . As a result, in the 
long wavelength limit A; — )■ 0, we have a result 1 + r{k — ?■ 0) = t{k — ?■ 0) = 0. From this 
condition, we have total reflection t = and r = — 1 in the limit A; — i- 0. We also obtain the 
phase shift of the amplitude reflection coefficient 9r = —iln (r/|r|) = vr. 
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Modes of (po and <^_i obey the Schrodinger-type equation. Generally, a particle obeying 
the Schrodinger equation shows total reflection against the potential barrier in the long- 
wavelength limit. The result of the excitation 0_i is consistent with this well-known result. 
The potential barrier makes the amplitude of the wave function vanish at zero momentum. 
However, the excitation 0o has a wave function with the same form as the condensate 
wave function, and its finite amplitude in the long-wavelength limit results in the finite 
transmission coefficient. This property is strongly related to the fact that the transverse 
spin wave is a Nambu-Goldstone mode. The barrier coupled only to the density does not lift 
the degeneracy of the ground state and it results in perfect transmission in the symmetric 
case 26!] . From another point of view, we comment that the wave function 0o behaves like an 
interacting Bose gas, where the potential effect disappears far from the potential barrier 
However, the wave function 0_i behaves like a free Bose gas, where the potential effect still 
remains far from it [2?^. These different properties lead to different results. 



VI. DEPENDENCE ON ci AND THE LINEAR ZEEMAN EFFECT 
A. Quadrupolar spin mode in the saturated magnetization state 

As seen in Fig. HI total reflection occurs in the long-wavelength limit in the quadrupolar 
spin mode in the saturated magnetization state; however, one finds total transmission at a 
small but finite momentum. We numerically investigate this high-transmission coefficient 
by changing the parameter Ci. 

We find that the tunneling properties of the quadrupolar spin mode can be categorized 
into two types. The first one is the excitation whose transmission coefficient touches unity 
(Fig. [T3|) . which we name type I. As seen in Fig. [13] (c), the phase of the amplitude reflection 
coefficient r jumps by vr at a certain momentum, which corresponds to the peak position 
of the transmission coefficients. This jump is a signature of type I, since the amplitude 
reflection coefficient crosses zero and changes its sign. When we focus on the momentum 
reaching T = 1 as a function of Ci/c+, the maximum k is found to exist from Fig. [13] 
(a). The momenta reaching T = 1 in Ci/c+ = —0.01 and —0.15 are less than that in the 
intermediate parameter Ci/c+ = —0.1. Type I can be seen in the small coupling constant 
regime of ci. The second type is the excitation whose transmission coefficients do not touch 
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unity (Fig, n^j) . which we name type II. As |ci| increases, the peak value of the transmission 
coefficient decreases from unity. 

In type I, the phase shifts of the amphtude transmission coefficient are 7r/2 in the long 
wavelength limit [Fig. [131(b)]. In type II, however, the phase shifts of the amplitude trans- 
mission coefficient are 37r/2 in the long wavelength limit [Fig. [T3] (b)]. A characteristic 
feature of the total reflection can be seen in the argument of the amplitude reflection co- 
efficient r. The argument of the amplitude reflection coefficient r reaches ivr in the long 
wavelength limit. Note that we can find type I and type II in the tunneling problem of a 
quantum particle obeying the Schrodinger equation, when we consider a mimic system of the 
present case (i.e., a system in the presence of attractive-repulsive-attractive potential) {28 1. 



B. Linear Zeeman effect 



In the saturated magnetization state, the uniform magnetic field associated with the linear 
Zeeman effect does not change a profile of the condensate wave function. As for excitations, 
it changes energy intervals between eigenmodes; however, it does not affect wave functions 
of excitations when one regards them as functions of the momentum and the position. As a 
result, the tunneling properties of excitations in the saturated magnetization state are not 
changed by the uniform magnetic field. 

In the unsaturated magnetization state, the density of each component is affected by the 
magnetic field associated with the linear Zeeman effect, as shown in Fig. [71 We thus have an 
issue regarding how the transmission coefficients are affected by the uniform magnetic field. 
From Fig. [151 one finds that the effects of the uniform magnetic field are quite small for the 
transmission coefficient, although the magnetization and the condensate wave function are 
easily changed by the magnetic field, as in Fig. [71 

A uniform linear Zeeman effect changes the spatial variation of the condensate wave 
function; however, the correspondence between the wave functions of excitations and the 
condensate wave functions still holds in the long- wavelength limit. As a result, the uniform 
magnetic field does not change the tunneling properties drastically, even in the unsaturated 
magnetization state. The tunneling properties depend on nonuniform effects, so that we 
expect that a nonuniform magnetic field changes the present results. We report this effect 
in a separate paper 
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From the viewpoint of conducting experiments, the Zeeman effect is useful to study the 
tunnehng problem of excitations in the spin-1 BEC. In the saturated magnetization state, 
the strong magnetic field can make spectra of three excitations (the Bogoliubov excitation, 
the transverse spin-wave mode, and the quadrupolar spin mode) separated well with an 
interval qubB. An excitation in a BEC can be stimulated by the Bragg pulse, an advantage 
of which is that one can produce an excitation with a particular momentum Kk and a 

r 

particular energy hoo using two laser beams with different wave vectors and frequencies [2£ 



3l|. If the energy levels of excitations are well separated in the spin-1 BEC, it is easy 



to stimulate each mode selectively making use of this resonance scattering. If one performs 
experiments on this tunneling problem of the excitation in the spin-1 BEC, it is interesting to 
investigate the saturated magnetization state. The Bogoliubov excitation and the transverse 
spin wave show total transmission; but the quadrupolar spin mode shows total reflection 
against the potential wall. Rich scattering properties of excitations could be seen in the 
saturated magnetization state of a spin-1 spinor BEC. 



VII. SUMMARY 



We determined the tunneling properties of excitations in saturated and unsaturated mag- 
netization states of a spin-1 BEC. In the saturated magnetization state, we found that 
anomalous total transmission in the low-momentum limit occurs in a Bogoliubov mode and 
a transverse spin mode. However, a quadrupolar spin mode undergoes total reflection in the 
same limit. In the unsaturated magnetization state, however, all excitations show anomalous 
total transmission in the low momentum limit. When the anomalous tunneling phenomenon 
occurs, the phase shift is always zero. 

In order to investigate what determines total transmission and reflection in the low mo- 
mentum limit, we focused on spectra of excitations and properties of wave functions. As 
for spectra of excitations, we have two points: the dispersion relation and the energy gap. 
As studied in Sec. IIIII for the saturated magnetization state, we found total transmission, 
irrespective of the dispersion relation. The Bogoliubov excitation, whose energy E is pro- 
portional to the momentum p in the low momentum regime, and the transverse spin wave, 
with E oc p"^, both show anomalous perfect tunneling in the low momentum limit. 

A key of anomalous tunneling was found in the energy gap. Indeed, we found that an 
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excitation without an inherent energy gap shows total transmission in the long-wavelength 
limit; an excitation with an inherent energy gap shows total reflection. We use the term 
"the inherent energy gap" as the energy gap in the absence of spatially uniform external 
fields. For example, it is not the energy gap due to the linear Zeeman shift induced by the 
uniform magnetic field. In the saturated magnetization state, the quadrupolar spin mode 
which shows total reflection has the inherent energy gap; however, the Bogoliubov excitation 
and the spin-wave mode which show total transmission do not have the inherent energy gap. 

Another key of anomalous tunneling was found in the properties of wave functions. We 
found that wave functions of excitations which show anomalous total transmission have the 
same forms as the condensate wave functions in the low- momentum limit. We note that the 
quadrupolar excitation showing total reflection has a wave function with zero amplitude in 
the long- wavelength limit, which is caused by the potential barrier. The spectra of excita- 
tions and the properties of wave functions suggest that the anomalous total transmission is 
strongly related to the existence of the zero modes in the absence of the uniform external 
field. 

Table [III summarizes the tunneling properties of excitations, their energy gap, and their 
wave functions in the low-momentum limit p — )> 0. Here, Ab=o denotes the inherent energy 
gap, and (pj denotes the wave function of excitations with j th component, Uj and Vj are 
the Bogoliubov- type wave functions given by (pj = Uj exp {—iEt/h) — v* exp {iEt/h) with an 
energy E, and $j denotes the condensate wave function with j th component. In addition, 
T and 6t, respectively, denote the transmission coefficient and the phase shift. 

We have some remaining issues. Here we studied the effects of the nonmagnetic potential 
barrier which couples only to the local density. It is unknown how a spin-dependent magnetic 
potential barrier affects the tunneling properties of excitations in spinor BECs. A spin- 
dependent magnetic potential barrier will drastically change the results of the nonmagnetic 
one. In the current carrying state of the scalar BEC, however, total transmission of the 
Bogoliubov excitation disappears at the critical current sj. In the current carrying state of 
the spinor BEC, it is also unknown what happens at the critical current for the tunneling 
problem. We report on these issues in a separate paper 
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TABLE II: Relations between the tunneling properties of excitations, their energy gap, and their 
wave functions. 

Saturated magnetization state (|m| = 1) 
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Appendix A: EXCITATIONS 
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Here we summarize excitations in the spin-1 BEG on the basis of earlier studies 
[19I . We will derive eigenenergies and eigenvectors from Eqs. ( ITT]) and ( IT2|) and discuss the 
properties of excitations. 
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1. Saturated magnetization state 



In the uniform system, equations for excitations around the saturated magnetization state 
($+1, $0, = (v^, 0, 0)T are given by (ESD, dMD, (EZD, and ([MD, with setting = 1, as 

/ 






V 2m 



2m 



V2 








(Al) 

(A2) 
(A3) 



where 5+1 = n+i + t>+i, and G+i = m+i — i;+i. 

For the spin component +1, the energy is given hy E = \J e{e + 2c+n), where e 
h'^k'^/{2m), and we have the following solutions: 



5+1 

Here, k and k are given by 



exp {±ikx) 



1 «+l 1 


h 1 











exp (±KX 




(A4) 



2m 



a/ (c+n)^ + i?^ — c+n 



2m 



v/(c+n)2 + + 



(A5) 
(A6) 



Coefficients a+i and /3+i, satisfying the normalization condition {a*^-^(3^i + q;+i/3^;^)/2 = 1 
(i.e., |m+iP — l^^+ip = 1), are given by 



a+i 



^/ (c+n)^ + i?^ + c+n 



v/(c+n)2 + E^-c^ 



n 



E 



(A7) 
(A8) 



For the spin component 0, the wave function 0o is a solution to a Schrodinger-type 
equation (\A2\i with the energy E = e + q^ibB. The solution in the uniform regime is given 



by 00 = exp (ii/cx) with hk = ^j2m{E — g/j^BB). 

However, (/)_i is also a solution to another Schrodinger-type equation (1A3I) with the energy 
E = e — 2cin + 2g^BB. The solution in the uniform regime is given by 0_i = exp {±ikx) 



with hk = >/2m [E + 2(cin — y/xe-B)]. 
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These fluctuations are related to the following modes 18|, ll9|]: </)+i is associated with an 
ordinary Bogohubov mode, which is composed of the density fluctuation 5n = <l'+i(0+i+0!5_i) 
and the phase fluctuation 69^i = — 0;^J/(2z$+i). (pQ is associated with the transverse 
spin wave mode (5M+ = -\/2yUB$+i0qjand 0_i is associated with the quadrupolar spin 



wave 



19| (the longitudinal spin wave 



18|) 5M2 



2. Unsaturated magnetization state 



In this state, the configuration of order parameters is given by Eq. (|H]). In the uniform 
regime Kxt(|a^| ~^ oo) = 0, we have the following equation: 
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S-i 
\G-iJ 
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2m 






2m 
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-D 
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2m 



-D 






\G-iJ 



where D± = c+n{l ± A) and D = c_a/1 — with A being gfisB /{cin). 
The energy is given by 

E{±) = ^/e{e + c+n±C), 



where e = fi?k'^/{2m) and C = a/ (c_n)2 + AcoCinPA"^. Solutions are given by 
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\g5V 



exp (±iA;(±)x) 















«(±) 













exp (±K(±)x) 



«{±) 

V "(±) / 



(A9) 



(AlO) 



(All) 



where the wavenumber k(^j--^ and the diverging and converging rate with the energy E 
are given by 



hk 



(±) 



ym" ^/{c+n±Cy + 4E^ - (c+n ± C) 



(±) 



m 



v/(c+n±C)2 + 4E2 + {c+n ± C) 



(A12) 
(A13) 
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Coefficients are obtained as 



«{±) 



2E{C±c^iiA) 



C[y(c+n ± C)2 + 4£;2 - (c+n ± C)] 
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2E{C T c+nA) 



(A14) 



\ V CW{c+n ± C)2 + 4^2 + (c+n ± C)] J 
which satisfy the normahzation condition Re[Q;(±)^(±) + q;(±)/5(±)] = 2, corresponding to 

E(kP-kn = 2. 
i=±i 

However, the equation of the spin component is given by 



E 



Uo 



' -— + cin 



2m 



cinA 
2m 



-V2 + cin 

The energy of this mode is obtained as 
Its eigenvectors can be written as 



Uq 
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— exp {±ikx) 



1 ao 1 
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f Mo 1 









— exp {±Kx) 



ao 



(A16) 



(A17) 



where the wavenumber k and the diverging and converging rate k are given by 

hk =^277?" V(Cin)2(l - A2) + £;2 - Ci77] , 



=^277?" V(ci77)2(l - A2) + £;2 + Ci77] . 

Here, ao and (3o are given by 

/ U{ciny{l-A^) + E^ + E \ 



(A18) 
(A19) 
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2E 



v/(ci77)2(l->12)+^2_^ 

2£; / 
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Eigenmodes with and respectively correspond to out-of-phase and in-phase 
modes between the spin components ±1, since signs of and S^^i (and also G^'* and 
G^_^i ) are opposite, and signs of and S^Si (and also G^^i and ) are the same. In 
Ref. 18|, these modes were termed the collective mode that couples with the number density 
and the longitudinal spin density. Specifically in the low energy limit, these excitations 
can be regarded as phase modes of condensates like the Bogoliubov mode. We thus term 
the mode with [-E(-)], "Bogoliubov and longitudinal spin excitation [out-of-phase (in- 
phase)]". The other mode 0o is associated with the transverse spin wave mode 6M^ = 
5M*_ = y2/iB($+i0o + $-i0o) bl- 



Appendix B: INTEGRABLE CASE 

As seen in Eq. fj42|) for the unsaturated magnetization state, the condition cq = ci makes 
the GP-type equation of the spin components ±1 decoupled. It also holds for the Bogoliubov- 



type equation 
leda et al. 32 



3]). The condition cq = ci corresponds to the integrable condition found by 
34| . We shall consider the special case cq = C2(= c). 
Using this condition, the GP-type equation is given by 



= ( -—V" + V,^t -CUT gf^BB $±1 + 2g$^i. 
2m 



We assume that the configuration is given by 



(<I>+1,<I>0,$- 



■ v/(l + A)n/2^+i (x) , 0, v/(l-AW20_i(x))- 



(Bl) 



(B2) 



where 0±i(|a;| — i- oo) = 1. Substituting it into the GP-type equation, we obtain the equation 







_ V' + Kxt(a;) + dn{l ± (7/iB5)(04i(x) - 1) 
2m 



>±i x 



(B3) 



This is the same form as the equation of the scalar BEG, and hence it is easy to solve. 
Far from the potential barrier, however, we have the equation of the excitation given by 



E. 



± 







-|^V2 + 2(c±^/XBi?) 
2m 







(B4) 



The energy is obtained as 



^± = \/e\e^2{in±g^^B)\ 



(B5) 
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and solutions are given by 

5±i 



G±i 



exp {±ikx) 



1 a±i 1 


h 1 











exp (±/tx) 




where 
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'2m 



cn ± gfiBBy + EI + (cn ± ^/xr^) 



and 



±1 



^(cn ± gjjBBy + El - {cn ± gfisB) 



E, 



We obtain the transmission coefficient by solving the following equations: 



-i—V + Kxt -CUT gi^BB + 2g$|i 
2m 

2m 




(B6) 



(B7) 
(B8) 



(B9) 
(BIO) 

(Bll) 
(B12) 



for a; — )■ — oo, (B13) 
forx-)-+oo. (B14) 

By comparing fIBll) and fIBlip for the limit -E — 0, we find that the wave function S±i 
in the long-wavelength limit has a solution corresponding to the condensate wave function 
^±1. We expect perfect transmission of both modes in the long wavelength limit. Figure [16] 
shows the results of the tunneling problems for the spin components ±1. We find the total 
transmission in the long- wavelength limit, and find that each phase shift goes to zero in 
this limit. In the integrable case cq = ci, we conclude that the fluctuations of different spin 
components are decoupled, and each excitation shows the total transmission in the long 
wavelength limit. 
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FIG. 8: Tunneling properties of out-of-phase and in-phase modes in the unsaturated magnetization 
state through the barrier 14xt(2;) = Vh^ix)- (al) The transmission coefficient T, (bl) the reflection 
coefficient R, and (cl) the phase shift of the out-of-phase mode for the case where the incoming 
mode is the out-of-phase mode. The transmission and reflection coefficients of the in-phase mode 
are not plotted. We found that these values can be regarded as zero, and one can find T + R = 1. 
(a2) The transmission coefficient T, (b2) the reflection coefficient -R, and (c2) the phase shift of the 
in-phase mode for the case where the incoming mode is the in-phase mode. We also do not plot 
the transmission and reflection coefficients of the out-of-phase mode in this case. We also found 
that these values are regarded as zero, and one can flnd T -|- i? = 1. We use the parameter of the 
magnetic field as qim-qB = 0.5cin, and the scattering lengths oq : a2 = 46 : 52, following Ref. [2]. 
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FIG. 9: (a) Transmission coefficients and (b) corresponding phase sliifts of tlie transverse spin- 
wave mode in the unsaturated magnetization state at B ^ 0. We use the same potential barrier 
and the same parameters as in Fig. [8j 
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FIG. 10: (Color online) Wave functions when the incoming mode is the out-of-phase mode in the 
unsaturated magnetization state. Normahzed functions £"±1 and respectively, are given by 

5±i = M^^^S±x and G±i = AA[t^G±i, where M^^^ = ^(1 ± A)l2^ECI\(C ± c+>l)(c+ + C)\. We 
use the (5-function potential barrier V(.yx{x) = V\j5{x) with Vb/Cp = 3con, and the magnetic field 
given by qhbB = O.bcin. Parameters oq : 02 = 46 : 52 are used, following Ref. [l^. 
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FIG. 11: (Color online) Wave functions when the incoming mode is the in-phase mode in the 
unsaturated magnetization state. NormaUzed functions S±\ and G+i, respectively, given by = 
a47^5±i and G±i = a4T^G±i, where AA^^^ = ^(1 ± A)l2^ECI\C ^ c+A)(c+ - C)]. We use the 
same parameters and potential barrier as in Fig. [TOl 
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FIG. 12: (Color online) Wave functions of the spin component in the unsaturated magnetization 
state. Normalized functions u and v are given by {u,v) = —^J A/n{uQ, vq). The potential barrier 
and parameters are the same as those in Fig. [TUl 
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FIG. 13: (a) Transmission coefficients, (b) corresponding phase shifts and (c) arguments of the 
amphtude reflection coefficients of the quadrupolar spin mode in the saturated magnetization state 
for the small coupling constant |ci|. For such a small parameter, the perfect tunneling occurs at a 
certain momentum, where the argument of the amplitude reflection coefficient jumps by vr. We call 
this type of property type I in this paper. We use the (5-function potential barrier Vext{x) = Vb5{x) 
with 14 /Cf = 5c+n. 
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FIG. 14: (a) Transmission coefficients, (b) corresponding phase shifts and (c) arguments of the 
amphtude reflection coefficients of the quadrupolar spin mode in the saturated magnetization state 
for the larger coupling constant |ci| than in Fig. [T3j In contrast to type I, perfect tunneling does 
not occur for a small momentum. We call this type of property type II in this paper. We used the 
same potential barrier as in Fig. [T3j 



43 




FIG. 15: Magnetic field dependence of the transmission coefficient in tlie unsaturated magnetiza- 
tion state through the (5- function potential barrier Vext{x) = Vh5{x) where Vh/^p = Scqti. (a), (b), 
and (c), respectively, show results of the out-of-phase mode, the in-phase mode, and the transverse 
spin mode. We use parameters oq : 02 = 46 : 52, following Ref. jl^. B is given in Fig. [71 
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FIG. 16: Tunneling properties of excitations in the integrable case (cq = ci). (al) Transmission 
coefficients and (a2) phase shifts of the spin component +1. (bl) Transmission coefficients and (b2) 
phase shifts of the spin component —1. We use the (5-function potential barrier Vext{x) = Vh5{x). 



We use parameters ao : 02 = 46 : 52, following Ref. 



19[. We also use g^sB = 0.5cn. 
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